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Abstract 

We introduce a method for generating explicit solutions to the non- 
linear Hodge-Frobenius equations, a large class of quasilinear variational 
equations. The method is specialized to construct solutions which are 
integrable in the Frobenius sense and applies to mixed elliptic-hyperbolic 
equations from geometric variational theory, shallow-water hydrodynam- 
ics, and singular optics. MSC2010: 35Q35, 35M10 
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1 Introduction 

Although there are many ad hoc arguments that work for particular equations, 
few methods are available for producing explicit solutions to entire classes of 
quasilinear partial differential systems. This is especially the case if the system 
changes from elliptic to hyperbolic type along a smooth hypersurface in its 
domain. We begin by considering the exact solvability of systems having the 
form 

V-(p(Q)w) =0, (1.1) 

Vxw = 0, (1.2) 

where w = (wi , . . . , w n ) is an unknown vector- valued function on a domain 
fi C E n ; Q = |w| 2 ; p : Dom p c K+ U {0} -> M+ U {0} is a prescribed, 
continuously differentiable function on its domain of definition. By eq. (jl.2p we 
mean, by perhaps a slight abuse of notation, the system of n x (n — 1) equations 

diWj — djWi = , for alH ^ j G {1, . . . n} , 
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equivalent to the requirement that all the 2x2 minors of the 2 x n matrix 



di d 2 ... d n 

Wl W 2 ... W n 



(1.3) 



vanish. Equations (11.11 II. 2[) are a vectorial special case (corresponding to the 
case of differential 1-forms in n dimensions) of the nonlinear Hodge equations 
introduced in |12j . The system is satisfied by a large number of physical and 
geometric models, including the irrotational steady flow of a compressible fluid; 
light near a caustic; shallow water hydrodynamics; ball lightning in the nonlinear 
conductivity model; and fields which can be represented as a non-parametric 
extremal surface in M 3 or M 3 . See Sec. 2.7 and Chs. 5 and 6 of [IT] , as well as 
references cited therein, for discussions of these applications. 

Equation (|1.1[) is a variational equation associated to the energy functional 
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E =o P^ ds * 1 ' 

z Jq Jo 

while equation (jl.2[) constrains the resulting field to be locally conservative. In 
the language of elementary differential equations, eq. (jl.2[) implies an exactness 
condition on contractible domains: if (|1.2[) is satisfied, then we can write 

w = VC (1.4) 

for some scalar function £. 

We develop the method of generalized streamlines, by which we derive the 
solution formula for the variational equation (jl.l[) . This formula holds globally 
on contractible domains and locally on more general domains. A less specialized 
version of the method, in which streamlines appear only locally, yields a global 
solution formula on non-contractible domains as well. Special extensions, by 
means of De Rham cohomology and Hodge decomposition theorems on mani- 
folds, can be found in [7], Sec. 5. 

Our method is qualitatively similar to classical techniques in which solutions 
to equations of fluid dynamics are constructed by following streamlines of the 
flow; see, e.g., Sec. 1.6 of [8]. But our method is more general, as it applies to 
equations which do not possess streamlines in the physical sense. 

Indeed, our method is general in the sense that it applies to the case of dif- 
ferential /c-forms presented in [12] . The extension to general fc-forms is straight- 
forward, as it is based only on derived consequences of the adjoint form of the 
Poincare Lemma, and every operation in the proof has an obvious analogy for 
forms of arbitrary order. This extension is illustrated in Appendix lAl 

Naturally, condition (|1.2j) - or (|4.2[) . in the case in which one is seeking 
solutions to eq. (|4.1|) . below - is not generally satisfied. So, continuing to 
follow the analogy of elementary differential equations, we look for an integrating 
factor. That is, we replace the constraint (|1.2[) by the Frobenius condition 

V x w = G x w, (1.5) 
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where G = G(x) is a vector function, or replace the constraint (|4.2I) by (|4.3[) . 
Again borrowing (and extending) the standard notation for the cross product 
of vector fields in R 3 , by eq. (jl.5|) we mean the system of n x (n — 1) equations 



diWj — djWi — GiWj — GjWi , for all i^j £ {1, 2, . . . , 



obtained by equating the 2x2 minors of the 2 x n matrix (11.31) to the corre- 
sponding minors of 



\ wi w 2 ... w n J 

As a consequence of the Frobenius theorem [3J , a solution to (ll.5[) in the context 
just described can always be written locally in the form 



where ry and £ are functions and G can always be chosen to be conservative: 



Thus one of the goals of the method is to construct explicit solutions to eq. (jl.ip 
which are completely integrable fields. By this we mean that they will possess 
an integrating factor, namely a function n satisfying (11.61) or equivalently, in this 
vectorial case, satisfying (11.51) . (|1.7[) . See [5J, Sees. 2.1, 2.2 for a discussion. 

Notice that the aim of our method is not to solve (jl.ll 11.51) for prescribed 
G, including the original case G = 0. Indeed, solving (jl.ll ll.5[) for prescribed 
G would entail imposing the condition that our explicit solution also solve an 
additional non-linear differential equation. However, a vector field G for which 
w satisfies (|1.5[> can be computed once the explicit solution to eq. (jl.ll) has been 
found. Moreover, as an artifact of the proof, we derive the auxiliary equation 
that would need to be satisfied in order to solve the system for prescribed G. 

Analogous arguments apply to the case in which one is considering eqs. 
I|4.ip . I|4.3p . below. The solutions to system (j4.ll) found by our formula are also 
completely integrable, as they ordinarily satisfy the Frobenius condition (j4.3l) 
for some vector function (11.71) . 

We apply the method of generalize streamlines to various typologies of sys- 
tems of equations, where the typology is determined by two different dimensional 
parameters, k (representing the degree of the form) and n (representing the di- 
mension of the domain). In all cases, the method decomposes into two steps: a 
first step yielding the normalized solutions, and a second step yielding the size 
of the vector solutions Q(x). A feature of the method is that the two steps are 
essentially unrelated. In particular, the first step is independent of the density 
function p and depends only on the typology of the system - that is, on the 
dimensional parameters k and n; the second step is totally independent of the 
system under consideration and is concerned almost exclusively with the den- 
sity function p, via the invertibility of the real- valued function of one variable 
£ = 4>{t) = tp (t), bearing in mind that the analysis of the singular sets will 
depend on the dimension of the domain. 




e -"w = vc, 



(1.6) 



G = Vt?. 



(1.7) 
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The system 11 II -5|) was introduced for the general case of differential forms 
of arbitrary order in [9] and studied in [10] ; it received an extensive treatment 
in |6j. We call that generalized form of (11.11 ll.5[) the nonlinear Hodge-Frobenius 
equations. Equation (jl.5|) represents the mildest natural weakening of the side 
condition on the variations in the conventional nonlinear Hodge equations; com- 
pare also with the conservative system studied in [J. The Frobenius condition 
(|1.5I) or (|4.3|) . and the special cases (11.2|) or (|4.2[) respectively, corresponding 
to the choice G = 0, generally have physical or geometric significance. For 
example, if the mass density p in eq. (|1.1[> satisfies 



then we obtain a model for steady, shallow, hydrodynamic flow under an appro- 
priate normalization; see, e.g., Sec. 10.12 of [14]. If (|1.2|) is satisfied, then the 
flow is irrotational; if ()1.5[) is satisfied, then the flow is vorticial (Sees. 15.21 and 
15.51) . If we choose in (jl.ll) the density 



then we obtain a model for extremal surfaces in 3-dimensional Minkowski space 
which can be expressed locally as the graph of a function £; see Sec. 1 of [13], 
Ch. 6 of [TT], and the references cited therein. Replacing (ll.2[) by (|1.5I) dilates 
the gradient vector of ( without changing the direction of that vector fSec !5.1[) . 
For other examples, see, Sees. 1-4 of [6]. 

1.1 Organization of the paper 

In Sec. [2] we develop the method of generalized streamlines in 2 dimensions for 
the system (Ti~Tl IT3]) . 

In Sec. [3]we extend the method to n dimensions, introducing the generalized 
stream flow generated by the stream matrix F. 

In Sec. [4] we apply the method to a related system of n(n — l)/2 equations 
in n dimensions. 

In Sec. [5] we present examples and applications to illustrate the cases con- 
sidered in Sees. [2HU 

In the Appendix we give a more general formulation of the method for dif- 
ferential forms of arbitrary degree on domains of arbitrary dimension. 

2 The method of generalized streamlines in two 
dimensions 

Denote by O a simply connected domain in R 2 . Applying the Poincare Lemma 
to eq. (jl.ip indicates that there is a scalar differentiable function / (x, y) such 
that 



p{Q) = 1 - | , 



(1.8) 



p(Q) = |i- QP 1/2 , Q^i, 



(1.9) 



P (Q)w = V±f, 



(2.1) 
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having denoted by Vj_/ the transverse gradient of / given by 



V x f = -(d y f)i + (d x f)j. 



In the language of fluid dynamics, eq. (|2.1j) asserts that there exists a stream 
function for the divergence-free vector field p(Q)w. Conversely, for any pre- 
scribed scalar function f (x,y) defined on a domain O, the field Vj_/(a;,y) is 
divergence- free. So (|2.1j) leads to solutions of (|1.1[> . even if f2 is not simply 
connected, in which case one may prescribe / with point singularities; cf. also 
Remark iv) at the end of this section. In seeking explicit solutions w of eq. (|l.l[) 
via (|2.1[) , one encounters two problems. The first is that eq. (|2.1j) may not 
be easily inverted to yield w in terms of /; the second is that p may vanish 
for values of Q attained on a subset 70 C £1. Both problems arise in actual 
applications. In what follows we illustrate how to proceed in order to invert 
eq. (|2.1|) and argue that in many cases the second problem does not occur, or 
is only apparent in the sense that it corresponds to the existence of removable 
singularities. 

Our initial observation is that, because p is a scalar function, eq. (|2.1[) yields 
for p 7^ the precise expression in terms of / for the normalized solution 

l w l |V/| 

The mass density p does not appear explicitly in eq. (|2.2I) . Nevertheless, this 
expression holds only in the limit (or not at all) on a singular set on which 
V/ = 0. Such a set would itself depend on p and /, as will be illustrated later in 
the present section and in the examples of Sec. [5j An immediate consequence of 
our observation is that, in order to find a complete expression for w, it suffices 
to find an equation satisfied by |w| alone in terms of /. To this end we take the 
squared norm on both sides of eq. (|2.ip and find that the real-valued function 
of a real variable 

4>{t) = tp 2 {t), t>0, (2.3) 

can be usefully inverted to yield t = Q = |w| 2 . 

We say that the function <f> restricted to an interval I\ or I 2 is invertible with 
inverse of type 1, or invertible with inverse of type 2, respectively, if 

|(Q)>0,Q6J 1 , or ^(Q) <0, QG/ ai . 

We denote by ipi, tp2 the corresponding inverses. The first choice corresponds to 
a region of il on which the system (jl.ll ll.5[) is elliptic; the second corresponds 
to a region of S7 on which the system 1 1 . 5|) is hyperbolic, as can be easily 
computed. In various models there are non-empty intervals 1\ and I2 on which 
<f> is invertible with individual inverses in the elliptic and hyperbolic regimes, 
although we may not know a priori the corresponding regions in the domain f2. 
Sometimes elliptic and hyperbolic solutions can be patched together along the 
sonic curve dividing the elliptic from the hyperbolic regime. We will illustrate 
some of these models in Sec. [5j 
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Because p is given and we have, by (|2.1[) . 

4>{Q) = Qp 2 (Q) = |V/| 2 , (2.4) 

the quantity Q is expressible in terms of the prescribed generalized streamline 
f (more precisely, in terms of the magnitude of its gradient) via the inverse, or 
inverses, of <f>. So a posteriori the elliptic and hyperbolic regimes correspond to 
regions of the plane where one uses inverses of type 1 or of type 2, respectively. 
Using (I2.ip . (|2.4I) to solve for w in terms of /, we obtain solutions to (11.11) 
expressed by 

V±/ ,„ 

w = p(V,(|v/P))< (2 - 5) 

where tjj may denote inverses on different monotonicity intervals. These can be 
defined on a subdomain 

Qf = {(x,y) 6fln Dom-yf : |V/| 2 (x, y) G Im 0} C Q , 

except possibly on a singular set S C 7 S ; here 

7s = {(x,y) €% : 4>' (V (|V/| 2 (x,y))) = or is undefined } , 

where one may switch between different inverse functions of <j). We observe that 
the set 7„ always contains the set 

10 = {(x,y)en f : p(^(|V/| 2 (x,y))) = and V (|V/| 2 (x, y)) ^ 0} , 

on which the solutions (|2.5[) may blow up or be undetermined. That is, 

7o C 7 S , 

and such an inclusion may be proper. This follows easily from the relation 

4>'{Q) = (Qp 2 (Q))' = p(Q) (p(Q) + 2Q P '(Q)) . 

Moreover, 70 C {(x,y) € Qf : |V/| 2 (x,y) = 0}. In fact, by §XZ§ - and its 
alternate form (12.61) . see below - the condition p (tp (|V/| 2 (x, y))) = implies 
|V/| 2 (x, y) = 0. This inclusion is proper when there are points (x, y) € Qf at 
which ip (|V/| 2 (x,y)) = and p (|V/| 2 (x, y))) ^ 0. Also note that, by the 
alternate form of (|2.5[) given by (|2.7I) below, it is clear that w is defined and 
equals zero at points where ip(\ V/| 2 (x, y)) = 0, even when p (ip (|V/| 2 (x, y))) = 
at these points. 

Although for a smooth prescribed function /, it may be possible to define a 
corresponding w in such a way to transition with continuity between different 
inverses ip — in particular, between elliptic and hyperbolic regimes - it will in 
general not be possible to do so with higher regularity. So, although w may 
be defined with some regularity on the set Qf — 70 , eq. (jl.ll) may not hold in 
the classical sense across the set 7,. For this reason, we regard 7^ as a possibly 



G 



singular set. In important applications, this set is a smooth curve - the sonic 
curve mentioned earlier. 

For analogous reasons, one also regards the set 700, defined as 

7oo = {(x, y) e ft : p (tp (|V/| 2 (x, y))) is undefined }, 

as a singular set for the equation. As / is prescribed, one can avoid altogether 
the sets 70 and 7,, and produce examples for which 70 = 7 S = and w is 
smooth. Nevertheless it is more interesting, mathematically and in terms of 
the applications, to produce examples for which 7 S and perhaps also 70 are 
non-empty. 

Strictly speaking, the set of points on which V/| 2 ) becomes unbounded 
does not belong to ftf, and the definition of w cannot be extended to these 
points, as is clear from the alternate expression for w given by (|2.7[) . 

For a given function /, one has ftf — ft - that is, the corresponding solutions 
w live on all of ft, except possibly on a singular set S C 7 S - if and only if 

H f = {\Vf\ 2 (x,y) : (x,y) € ft} Clmc/). 

In the particular case in which 

Hf C Im <j) , 

the full range of <f> is not available and one may be able to invert only in the 
elliptic regime, or only in the hyperbolic regime, for that particular choice of 
/. There are also choices of / for which ftf — 0. The example of the density 
function p = k/^/Q illustrates this possibility when |V/| 7^ k. In fact, for this 
density function one has 

-4w = Vx/, 



which admits solutions if and only if |V/|(a;, y) — k V(x,y) G ft. Notice that 
only the normalized w/y/Q is specified by this equation, while the function 
Q(x,y) can never be determined. In fact in this example 4>{Q) = k 2 cannot 
be inverted and one always has ~ 70 C 7 S = ft. In the case \Vf\(x,y) = k 
\/(x,y) G ft, one has ft f — ft and the normalized solution corresponding to 
/ is defined everywhere. Multiplying by the arbitrary function y / 'Q(x, y), one 
obtains a family of solutions. This dramatically illustrates the importance of 
the range of |V/| 2 in determining solutions to eq. (Il.lj) . 

The simple example of the density function p = illustrates instead a 

case where 700 ^ 0, and yet one can still carry out our method and obtain 
solutions 

w = Vi/ |V/|, 

defined on the domain of V/ and smooth for smooth functions /. In this 
example, <f> = ^/Q is everywhere defined and increasing. Nonetheless, although 
the solutions defined above may be everywhere smooth, they are technically 
solutions to only on {{x, y) € ft : \7f(x,y) ^ 0}. 
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Remarks. 

i) In (|2.5p we have divided by p - which may vanish - but, as mentioned 
previously, this problem is only apparent in many cases. In fact, suppose that 
p(Q) = for some finite value(s) of Q, say for example p{Q\) = 0. This gives 
a posteriori that |w| = y/Qi is bounded at the point (s) pi £ fi at which the 
value Qi is attained, despite the vanishing of p at p\. It is possible, however, 
that in some cases w may be not well defined at p\ ; cf. Remark ii below. That 
is the case for example in some models in which w develops point singularities. 
It is also possible that p will vanish at points of fi for which Q becomes infinite. 
That is the case, for example, in the variational model considered in Sec. 15.11 
when one chooses to invert in the hyperbolic regime; see eq. (|5.5[) for |V/| — > 1. 

ii) Rewriting (|2.4p as 

|V/| 2 = 0OMIV/I 2 )) = p 2 W|v/| 2 )) VX|V/| 2 ) , (2.6) 

we obtain an alternate expression for w which may be used in place of (|2.5j) for 
finding explicit solutions to (jl.ll) . namely, 

Vi/ 



w = ±^0M|V/| 2 ). (2.7) 

This illustrates what may happen at points of 70 : even though Q may remain 
bounded at these points, in various examples and applications the normalized 
vector function w/|w| may not be defined (even as a limit) at these points,. 

Formally, the two formulas for w, (|2.5p and (12.71) . coincide if the minus sign 
in (|2.7p is assigned to regions where p (?A(|V/| 2 )) is negative (generally excluded 
in the applications). On the other hand, the minus sign in (|2.7p may be avoided 
altogether, as 

|V/|(z,y) = |V(-/)|(z,y) ; 



the choices of either — / or / in (|2.7p are equally legitimate and yield the same 
range: E_/ = 

Hi) The vector function w given by eq. (|2.5p or (|2.7|) satisfies (jl.ljl but may 
not satisfy (|1.2p . However, for any given G, including the case G = 0, eq. (jl.5p 
will be satisfied with the choice 

G Gi wwm ' ( } 

whenever / satisfies the equation 

A/ + Gi-V/ = (2.9) 



for some continuously diffcrcntiablc vector function Gi. Solving (|2.8p for Gi and 
substituting the result into (|2.9p , we find that the latter is a nonlinear equation 
for /. In any case, our aim is not to solve eq. (|2.9p for / in terms of a prescribed 
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Gi, but to obtain completely integrable solutions by the method of generalized 
streamlines. 

Indeed, for the reasons explained in Sec. 1, we require only that condition 
(|1.5[) be satisfied for some G once w has been determined by formula (|2.5[) . 
This can always be achieved, except on the (possibly empty) set 

7G = {(x, y) G fl f C M 2 : Vf(x, y)=0; Af(x, y) ? } , (2.10) 

by choosing for example, 

Gl = - a/ (w)' (2 - n) 

Because adding to Gi a vector function H satisfying H • V/ = has no effect 
on eq. (|2.8[) . there will be, in general, infinitely many G satisfying eq. (|1.5[) . 
We know from the theory [3] that one can specify G to be a conservative vector 
field. 

iv) The solution formula to (jl.ll) found by the method of generalized stream- 
line can be extended to domains f2 which are not simply connected, by replacing 
the transverse gradients V j_/ in (|2.1[) by arbitrary divergence-free vector fields 
a. Once this substitution is made, the remainder of the procedure does not 
change, yielding 

a 

w = 



pWH 2 )) ■ 

A calculation analogous to the one done in Remark Hi) shows that the Frobenius 
condition (|1.5[) for w holds with G satisfying (|2.8[) and Gi given by 



Gi = -(9 2 Q!i - dia 2 ) ( ) = ~ V -L ' a ( t-t 2 



with Vj_ = (— 82,81). 



3 The method of generalized streamlines in n 
dimensions 

The method described in Sec. [5] for the 2-dimensional case can be extended to 
vectors in n dimensions as follows. 

The Poincare Lemma applied to eq. (jl.ll) on contractible domains 57 C W 1 
gives the existence of an n X n real matrix- valued function F = (fij), with F 
skew symmetric (fa — and fij — —fj%, Vi,j € {1, . . . ,n}), such that 

p(Q)w = V ± F, (3.1) 

where by analogy with Sec. [2] we define the transverse gradient of a skew sym- 
metric matrix F, denoted V ±F, to be the n-vector 



V_lF 




nj I L n ■ 
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(Notice that for n = 2 this corresponds to the description in Sec. Rafter making 
the substitution / = — f 12 = /21.) 

Extending the language of fluid dynamics, eq. (|3.1[) asserts that there exists 
a stream matrix F for the divergence-free vector field p(Q)w. On domains Q 
which are not simply connected, eq. (|3.1|i still leads to solutions of (|1.1[> . More 
generally, one can replace Vj_-F by arbitrary sufficiently smooth divergence-free 
n-dimensional vector fields a on the right-hand side of (|3.1I) . 

In order to find explicit solutions w to eq. (jl.lj) in terms of given matrices 
F, we continue to proceed by analogy with Sec. [5J 

It is again possible to express the normalized solution to (jl.ip in terms of 
the stream matrix by the formula 

|w| |Vj_F| 

in which the particular choice of mass density p does not appear, although the 
singular set for this formula may depend on p. 

Because p is given, eq. (|3.1I) yields as in the two-dimensional case an explicit 
relation between the two scalar quantities |Vj_F| 2 and Q = |w| 2 , namely 

<P{Q) = Q P 2 {Q) = \V±F\\ (3.3) 

enabling us to express the quantity Q in terms of the prescribed stream matrix 
F via the inverse(s) of (j>. 

As in the two-dimensional case, the elliptic and hyperbolic regimes for the 
system of eqs. (|1.1[) . (|1.5|) correspond to the regions of ft C R" on which one 
uses inverses of type 1 or of type 2. Denoting simply by ip the various inverses 
of 4> restricted to its monotonicity intervals, one obtains solutions to (jl.ll) of the 
form 

V 1 F 

w = — (3 4) 

These are defined on subdomains 

n F = {x= (x 1 ,. ..,x n ) e nn Dom v±F : |V_lF| 2 (x) G Im<j>} C Q,, 

except possibly on a set S C 7 S with 

7s = {x G fl f : cf>' (ip (\\7 ± F\ 2 (x 7 y))) = or is undefined } . 

We observe once more that the set 

7o = {xeO f : p(^(|V ± F| 2 (x))) = and VXI V ± F| 2 (x)) ± 0} , 

where w may blow up or be undefined, is contained in j s . The latter is to be 
regarded as a possibly singular set. Even though, for a smooth prescribed stream 
matrix F it may still be possible to define w as a vector function with continuity 
on the set dividing the elliptic from the hyperbolic regime, eq. (11 . 1 1) will not 
hold on this set. In this context, under suitable regularity assumptions, the sets 
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70 and 7 S will be smooth (n — l)-dimcnsional hypersurfaces. As F is prescribed, 
one can keep away from these singular sets and produce various examples for 
which 70 = 7 S = and w is smooth. Nevertheless, it is again more interesting, 
intrinsically as well as in terms of possible applications, to produce examples 
for which 7 S , or 70 and 7 S , are non-empty. Again, the set of points for which 
tp(\W F\ 2 )(x,y) becomes unbounded is excluded from £1/ by definition, and w 
cannot be extended to these points, as is clear from the alternate expression for 
w given by 

w = ±^!^(|V^| 2 ), (3.5) 

an n-dimensional analogue of (|2.7p . 

For a given stream matrix F, one has £If = H - that is, the solutions w 
given above live on all of f2, except possibly on the sets 70, on which w may be 
undefined, and 7 S , on which w is not technically a solution to eq. (jl.lj) - if and 
only if 

E F = {|V±F| 2 (x,y) : (x,y) e ft} C Im<£. 

In various models there are non-empty intervals I\ , I2 on which <j> is invertible 
with individual inverses on the elliptic and hyperbolic regimes. Some of these 
models will be reviewed in Sec. 5. 

The analysis of the Frobenius condition is similar to the analysis carried out 
in Sec. O The Frobenius condition is obtained with 

G -- Gl -^mWiWT' (3 - 6) 

where Gi satisfies the system of n(n — l)/2 equations 

V x Vj_F + Gi x V ± F = . (3.7) 

On domains f2 which are not simply connected, in the more general context the 
Frobenius condition holds with Vj^F in eqs. (|3.6p . (|3.7[) replaced by arbitrary 
sufficiently smooth divergence-free n-dimensional vector fields a. 



4 The method of generalized streamlines for sys- 
tems 

The method of generalized streamlines also applies to the problem in which eq. 
p.ip is replaced by the system of equations 

di(pwj) — dj(pwi) = , for all !/je{l,...n}, (4.1) 

for a vector function w(x) = (wi, . . . , «j n )(x) on a domain £1 C M™, for pre- 
scribed p = p(|w| 2 ). Equation (|4.ip prescribes the vanishing of the 2x2 minors 
of the matrix 

/ di d 2 ... d n \ 
V PWi pw 2 ■ ■ ■ pw n J 
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In this case, eq. (I1.2j) and its more general version (|1.5[) are replaced by 



Vw = VdjU>j = 0, (4.2) 



3=1 



and 



w = £ a i w J = G ' w = £ G 3 W J » ( 4 - 3 ) 



respectively, where G(x) = (Gi, . . . , G„)(x) is a vector function. 

The method of generalized streamlines applied to this problem yields the 
solution formula for classical solutions on contractible domains 

V/ 

w = p(V(|v/P))' (44) 

or the alternate formula 



|V/| 



VW/F5). (4-5) 



in which / : fi C M" — >• K is any sufficiently smooth function and ^ denotes 
the inverse(s) of the function <f> defined by (|2.3I) . Vector fields w expressed by 
(|4.4p or (|4.5j) are defined on the set fl t except possibly on a singular set S C j s . 
Often they can be defined with some regularity on ftf\jo, although (|4.1[) is not 
technically satisfied on 7 S ; cf. Sec. [2] 

In order to obtain solution formulas on non-contractible domains f2, one can 
replace V/ in (14.4[) and in (|4.5p by arbitrary sufficiently smooth vector fields a 
defined on fl, satisfying djCti = diCtj V 

The Frobenius condition eq. (|4.3[) can be obtained for vector functions G 
satisfying the equation 

1 V(pOH|V/| 2 ))) , , 

A/ -2 V/ - P()(|WP)) - G - V/ ' (4 ' 6) 



V • a - - a • Vttt = G • a (4.7) 



or 

v(pWH 2 ) )) 
2" pOMM 2 )) 

if V/ is replaced by a as specified above on a non-contractible domain O. 

This section extends the results of Sees. [2 [3l to the system (|4. 1|) . The 
results in Appendix A can be applied to extend the method to more systems. 
For an application of the contents of this section, see Sec. 15.41 
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5 Examples 



5.1 A quasilinear elliptic-hyperbolic variational problem 

We first illustrate the method on a singular, quasilinear system of variational 
equations corresponding to II .5[) with p given by eq. (II. 9[) . The resulting 
elliptic-hyperbolic system provides a model for certain nonparametric extremal 
surfaces, as observed in Sec. 1. 

Suppose that a surface £ is locally the graph of a function £ with w = V£ 
satisfying eqs. (jl.ll ll.2|) with p given by (jl.9[> . Replacing condition (jl.2[) with 
condition (| 1 . 5|) corresponds to multiplying, at each point of TE, the length of 
the gradient vector V£ by a conformal factor exp[ry]. This corresponds to a 
conformal dilation in the noneuclidean part of the local metric tensor 

a- — T ( ^ x 

where / is the 2x2 identity matrix corresponding to the trivial (euclidean) 
metric. In this example, the function </> evaluated at Q is 

and 

dHQ) _ , i 

where the plus sign corresponds to the interval Q < 1 and the minus sign to the 
interval Q > 1, respectively. For Q < 1, which in the case G = corresponds 
geometrically to a maximal space-like hypersurface, the system of equations 
11.51) is elliptic. In that case, the function ip B = [cj)\ ]~ : [0, oo) — > [0, 1) 
is given by 

V's:^^. (5.2) 

For Q > 1, which in the case G = corresponds to a maximal time-like hy- 
persurface, the equations (jl.Hll.5p are hyperbolic. In that case the orientation- 
reversing function 0| is again invertible, with inverse given by the orientation- 
reversing function tpt — ) ■ (1>°°) (1)°°); 

i'f-C^^j- (5.3) 

For any given function /, one may always use ip = ip s as defined in (|5.2p . 
substitute into eq. (|2.5p . and obtain a solution to (II. ip : 

w s = - V±/ . (5.4) 
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This solution is at least as regular as V/ and remains bounded - in fact, Q 
remains in the interval [0, 1) . We observe that, although such w s is smooth and 
bounded on R 2 , p becomes singular over the set 

7oo = {(x ,yo) €tt 2 : lim |V/| 2 = oo| 

corresponding to Q — > 1. Of course this set will be empty provided / is prescribed 
so that V/ is bounded on bounded domains. If 7^ 0, the vector function w 
may still be defined as a limit. For these solutions 70 = and 7 S = 7^; cf. Sec. 

m 

If |V/(a;, y)\ £ [0, 1), the use of (15. 2[) to invert is the only possible choice. 
However, on regions of R 2 on which \Vf(x, y)\ £ (1, 00) one may use either (|5.2p 
or (|5.3p to obtain solutions to When the latter choice is made, one obtains 

w ' = 7mht (5 - 5) 

As in the previous case, although such vector functions W( remain bounded 
and may be well-defined on the set 700 (because Q — >• 1 for |V/| — > 00), eq. 
(jl.ip is not satisfied at 700, as p would become singular. In distinction to the 
previous case, w t becomes unbounded if |V/| = 1 on some subset of Q. In this 
case fi/ C is a proper inclusion as the set on which Q — ip(\^7f\ 2 ) becomes 
unbounded is technically excluded from the definition of O/; cf. Sec. [5] 

The choice to invert <j> using (|5.2[) throughout corresponds to a choice to 
solve the system (jl.ll 11.51) in the elliptic regime. A solution is obtained except 
on the set 700 . Alternatively, one may choose to use (|5.3[) in one or more regions 
on which \V f{x,y)\ £ (l,oo). It is easy to see that the 1-forms w s and Wj may 
only be patched together to make a continuous vector function w along curves 
in 7oo, as Q — > 1 for (x,y) approaching 700. The resulting vector function w 
would solve the elliptic- hyperbolic system 1 1 . 5|) except at points of 700 . In 
order to avoid singularities for w, one should choose / so that |V/| is bounded 
away from 1 on bounded regions on which we choose w = w t . 

5.1.1 The patching 

As a simple example to illustrate the patching along 700 of a solution to system 
(|1.1[) of type w s with a solution of type w t , with p given by (|1.9p . we consider 
the radial function f(x,y) = ±log|r — 1|. Here r = x 1 + y 2 , where the plus 
sign holds for r < 1 and the minus sign holds for r > 1. This function and its 
transverse gradient, 

V ^ = ^TT (--<-)' 

\r — 1| \ r r v 

are defined except on the set {r = 1} C M 2 and at the origin. The corresponding 
solutions Ws and Wj, expressed by (|5.4p . and (|5.5p respectively, are 

WS; * ~ Jl±(r-1) 2 V r' r) ' 
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The solution w s is denned on R 2 \{(0,0)}, although / was not defined on the 
set {r = 1}, while w t is defined on {(x, j/)eR 2 : r £ (0, 2)}. The two solutions 
can be considered separately, or one can obtain additional solutions by patching 
along the circle 

loo = {(x,y) eR 2 : r=l}, 

as on this circle w s = w ( = {—y/r,x/r) . An elementary calculation shows that 
this can be done with continuity of the first partial derivatives. For example, 
one can define w = w t on the punctured disk {(x, y) el 2 : r € (0, 1)}, and 
w = w s on the external region {(a;, y) £ M 2 : r € (1, oo)}. The resulting vector 
function w satisfies w e C 1 (R 2 \{(0,0)}) . Technically w solves (jl.ip . with p 
given by (|1.9[) only on the set R 2 \{r = 1}. A direct calculation employing 
eqs. (|2.8[) and (12.9[) shows that w satisfies the Frobenius condition (jl.5l) with 
G conservative. For such G and the given choice of p, the vector function w so 
constructed is a solution to the elliptic- hyperbolic system (11.11 II. 5|) . 

5.2 Shallow flow with vorticity 

Recall from Sec. 1 that if we choose p as in eq. (|1.8|) . then we obtain from eqs. 
([1.11 ll.2[) the continuity equation for shallow, steady, irrotational hydrodynamic 
flow under a convenient normalization. In this case w is the flow velocity. The 
flow potential £ exists locally by condition (|1.2I) and satisfies 

p(Q)d x (=-d y f, p(Q)d y ( = d x f, (5.6) 

where / (x, y) is the stream function. That is, in this application the function / 
of SecJ^Jhas a physical interpretation as an actual stream function; but if (11.21) is 
replaced by (| 1 . 5|) , the local flow potential £ exists modulo an integrating factor 
and the methods of, e.g., [1], Sees. 8 and 9, do not apply in an obvious way. In 
that case eq. (jl.ip is transformed locally into the system 

p(Q) wi = ~d y f , p(Q) w 2 = d x f. 

Multiplying this equation by e~ v and imposing the Frobenius condition (|1.5p 
with G = V17, eqs. (I5.6P are replaced by 

p(Q)d x C = -e~ r 'd y f, p(Q)d y ( = e-vd x f. (5.7) 

Equations (|5.7p allow for an interpretation of / and C as A-harmonic functions, 
that is, as functions satisfying the generalized Cauchy-Riemann equations 

A(x,|VC| 2 )-V ± /. (5.8) 

In this case, A(x, |VC| 2 ) = p (e 2 ^ \V (\ 2 ) VC ■ 

Proceeding as in Sec. [21 we find that the function ef> is strictly increasing 
on the regime of tranquil flow, on which the equations ([1.11 1 1 . 5j) are an elliptic 
system, and strictly decreasing on the regime of shooting flow, on which those 
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equations are a hyperbolic system. More precisely, in this example the function 

<KQ)b 

Q N 2 



<t>(Q) 



l 



(5.9) 



and its derivative is 



= (1-?)(J-^) 



which is strictly positive for Q e [0, |) U (2,+oo) and strictly negative for 
Q E (§7 2). However, when the squared speed Q exceeds 2, the mass density p 
becomes negative and the physical model no longer applies. Because 



and 



0(2) = 0(0) = 0, 



we have the following correspondences of explicit ranges on which the function 
4>{Q) can be inverted: 



(5.10) 
(5.11) 



If we are willing to consider negative values of p, then we also have a third case: 

(2,oo) -> (0,oo). 

The correspondence (I5.1ip is orientation-reversing and corresponds to the hyper- 
bolic regime for the system (jl.ll) . (I1.5p . The explicit expression for the inverse 
of on each of the three ranges above is given by the cubic formula. Denote by 
ipl, ip2, V>3 the three inverses of (f> corresponding to Q € [0, |), Q € (| ,2), and 
Q G (2,oo), respectively. When 

|V/Pe(o, (§)*), 

the function 0> is not uniquely determined, but may be chosen to be tpj for 
j = 1,2, or 3. If 

|v/P>g) S , 
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then the only possible choice is ip — -03- Obviously, when making these choices 
care must be taken to obtain reasonably smooth solutions w which have a 
plausible interpretation in the context of the chosen model. Note also that the 
above intervals depend on the constant in Bernoulli's law, which affects the 
constants in (|1.8[) : see the discussion of eq. (2.43) in [TT] . 
We will return to this example in Sec. 15.51 

5.3 Light near a caustic 

If we choose 




(5.12) 



then eqs. II. 2[) arise in an elliptic-hyperbolic model for light in the neigh- 
borhood of a caustic with the constant r 2 representing the index of refraction; 
see Sec. 1 of [5], Sec. 3 of [2J, or Ch. 5 of [TT]. The function <f) is invertible on 
the shadow region of attenuated waves, which corresponds to the elliptic regime 
of the system - that is, to the interval I\ = [r 2 ,oo). The function <f> is again 
invertible on the illuminated region of propagating waves, corresponding to the 
hyperbolic regime of the system - that is, to the interval I2 = [0, r 2 ]. The bound- 
ary between the regions of darkness and light is the caustic, on which the type 
of the system (jl.ll II. 2[) degenerates. As in the preceding cases, replacement of 
(jl.2l) with (jl.5l) introduces a multiplicative conformal factor in eq. (| 1 .4[) which 
complicates standard techniques, such as quasiconformal mappings, which are 
based on the existence of a field potential; cf. [2], Sec. 3. 

Applying the arguments of Sec. [H we obtain the functions 

i> d : [0,oo) -> [r 2 ,co), £^£ + r 2 , 

ipi : [0, t 2 ] — > [0, t 2 ] , £ — > t 2 — £ (orientation reversing) , 

yielding 

w=(V^|W)-J^, (5-13) 

where the plus sign corresponds to the choice ip = ipd, and the minus sign, to 
the choice ip = ipi, in (|2.5[) . We conclude from this formula that the condition 
I V/ 1 2 7^ is not needed for various choices of /, although relaxing this condition 
may lead to singularities (in particular, to point singularities when |V/| = at 
an isolated point). Again, for any given / one can use ipd throughout and obtain 
a solution w (corresponding to the plus sign in () 5 . 1 3 [) ^) . defined everywhere 
except possibly at points at which |V/| = 0; one may also decide to use tpi in 
one or more regions of the plane at which |V/| 2 < r 2 . The two types of solutions 
can be patched together only along curves 7 on which |V/| =0. If / is chosen 
to satisfy (|1.5I) for some conservative vector field G, then one obtains in this 
way a solution to the elliptic-hyperbolic system (11.11 11.51) . 



17 



5.4 The Born Infeld model 



In this example we consider the system of equations (|4.1[) in 3 dimensions, with 
p(Q) as in (|1.9p . As computed in Sec. 2] for the general case of n dimensions, 
the solutions w can be expressed locally by 

V/ 

w pWIv/I 2 )) ' 

and we will find below explicit forms for the expression p(-0(|V/| 2 )). 

For a physical interpretation of this system with the additional equation 
(|4.2p with n = 3, see Sec. 1 of [13] ■ For the effect on the model of using the 
above equation as a replacement for the more general equation (|4.3I) with n = 3, 
see Sec. 4 of [5J. 

This model and the one studied in Section 15.11 although very different, use 
the same type of density function p. The two models are dimensionally different: 
the one studied in Section \5. II is really a model for 1-forms in 2 dimensions, while 
the one studied here is a model for 2-forms in 3 dimensions. As a consequence, 
the equations considered are different for the two models. Nonetheless, as the 
density function p and the function </> are the same in the two models, much 
of the analysis carries over to the present model. So, as already computed 
in Section 4>(Q) = Q/\l — Q\, with Q / 1, admitting the inverses ipi = 
[^l[o,i)] _1 : [0,oo) ^ [0, 1), and i/> 2 = : (1, oo) ->■ (1, oo), given by 



We obtain the corresponding solutions 



£±1 



Wi-2 = — p= { (5.14) 

The solutions wi are bounded and defined for every choice of generalized stream 
function / (of class C 2 ), while the solutions W2 are defined for generalized stream 
functions / (of class C 2 ) such that |V/| 2 > 1 and become unbounded if |V/| 2 = 
1 at points of the domain Q. In principle, one could also prescribe a stream 
function / such that |V/| 2 — > oo when approaching a surface, say 7oo, contained 
in il. Solutions of type Wi and of type W2 could then be glued together along 
700. Nevertheless, the equations (|4.1j) with p prescribed as in (|1.9[) would no 
longer be satisfied, as p would blow up at the points of . 

In order to find vector functions G so that eq. (14.31) will be satisfied as well, 
we impose the condition 

V A/ 1 V(|V/| 2 )-V/ _ G-V/ 

VlV/Pil 2 (|V/| 2 ±1)^ vWF±i' 



Multiplying through by \/|V/| 2 ± 1, we obtain an equation for G of the form 

1 V(|V/| 2 ) 

U-tWfirtT -V/ = G.V/. (5.15) 
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From this we see that the relation (|4.3p can be satisfied with 



2 v1W±l l y /l 2 ' 

defined on the domain of the solution, with the possible exception of the set 

7G ee {(x, y, z) e Dom w C K 3 : Vf(x, y, z) = ; A/(x, y, z) ^ } , 

in the case of solutions of type wi. Note that if G is a solution to eq. (|5.15l) 
for given /, then the vector functions G + H with H • V/ ee are also solutions 
to eq. (|5.15|) . One needs to use this degree of freedom to find a conservative G 
satisfying (|4.3|) . 

In particular, choosing / to be the fundamental solution of Laplace's equa- 
tion in dimension 3 

/ = f(r) = -, r= y/x 2 +y 2 + z 2 
r 



yields the solutions 



^ 2 = — r%/1±r4 ' ( 5 - 16 ) 



which satisfy the Frobenius condition (|4. 3[) with 

G = 2^±P±^+H. 

r 4 Vl ± r 4 

The vector field G can be made conservative by a particular choice of H satis- 
fying H ■ V/ ee 0. 

As observed earlier for the general solutions W2, the solution (|5.16[) corre- 
sponding to the minus sign becomes unbounded on the set 

{(x,y,z)eR 3 : |V/| 2 (z, y, z) ee - = 1} , 

that is, on the sphere r = 1. This can often happen when inverting in the 
hyperbolic regime. In fact, in that case the map cf> is orientation-reversing and 
Q can become unbounded at points on which the prescribed |V/| is bounded, 
and thus on a bounded domain. 



5.5 Special classes of streamlines 

By prescribing certain symmetries in the stream function / {x, y) , one obtains 
explicit examples of solutions w to (jl.ip . For example, assume that the function 
/ can be expressed as f(x, y) — f(t), where t — g {x, y) is a given continuously 
differentiable function of the given domain f2. We obtain 

d x f = f(t)d x g, d y f = f'(t)d y g, 
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and 



For any choice of p, eq. 



IW| 2 = [f'(t)\ |Vg| 2 . 
yields a solution w of p. II) satisfying 



w 

Iwl 



± 



A special case is the example of radial / - that is, t = g (x, y) = x 2 
for any choice of p we obtain the simple formula 



y 



2 . Then 



w 

Iwl 



±- 



-yi + xj 



(5.17) 



In this case formula (12.41) becomes 



and (|2~7l) yields 



w = ± 



HQ) = |v/| 2 = uf (t) 

yi + xj 



Vt 



In applications to the shallow- water model corresponding to (|1.8[) . one may 
choose the function / so that |V/| 2 (t) is the function cf>(t) modulo a constant 
dilation of the domain. In that case, ip (|V/| 2 ) is essentially the identity on the 
domain. We illustrate this application with a numerical example. Choosing 



f(x,y) = [t 



t 2 



1 



where R is a given positive constant, we have 



|v/| 2 W 



i 



t 



t 



2RJ R 

We obtain, as expected, a solution having the familiar form 



^(^(l"^) 2 )/2 



I 1 2r) 1 + Jn ( 1 2fl) 3 _ -yi + xj 
1 - JL _ 

x 2R 



R 



where we have applied the various inverse functions ip described in Sec l5.2l on 
the various subdomains B 2 r/3, B 2 r\B 2 R/3, and R 2 \B 2 r- In the second to last 
equality we have used the identity 



i> p 



t \ t 

RJ R 



if) O (f> 



t 

R'' 



t = x 2 +y 2 , 
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which is satisfied by all the inverses. 

The vector w so defined satisfies with p given by (| 1 .8[) - that is, with 



and (IT3t with 



t 

2~R' 



G = — = j H = j 

on R 2 \{(0, 0)}, and is globally smooth. The system (jl.ll II. 5p is elliptic on 
^2_r/3\{(0, 0)}, hyperbolic on B2 r \B 2 R/3, and then again elliptic on R 2 \i?2_R 
(although that region would correspond to negative p). 

More examples can be found by applying the same procedure to a more 
general function t = g(x,y), such that |Vp| would be a function only of t. 



A The general case of &;-forms in n dimensions 

The extension of the method of generalized streamlines to differential forms of 
arbitrary order in arbitrary dimension is immediate. In fact, the method can 
be described naturally in the unifying language of differential forms; see [7] for 
a development in this direction. In this more general context the nonlinear 
Hodge-Frobenius equations appear in the form 

5{p(Q)u) = 0, (A.l) 

duo = T A duj, (A.2) 

where we solve for uj, a differential fc-form on a domain Slcl";d: A k (Q) ->■ 
A fe+1 (f2) is the exterior derivative with formal adjoint S : A k (fl) — > A fe_1 (f2); 
Q = |u;| 2 ; p : Dom p C K+ U {0} -> M + U {0} is a prescribed, continuously 
differentiable function on its domain of definition; and T is a 1-form. 

Condition (IA.2I) guarantees complete integrability in the cases k = 1 and 
k = n— 1. If r can be made exact, say T = drj 1 uj is said to be gradient-recursive 
and can be written as 

u = e"rfC 

for a (k — l)-form £; cf. Sec. 2.2 in [B]. 

If the domain fl is contractible, then applying the Poincare Lemma as in 
Sees. [5J [31 or 21 we obtain the solution formula for fc-forms uj satisfying (|A.ip 
in terms of generalized (n — k — l)-stream forms /: 

• ^ (A.3) 



where * : A* — > A" k is the Hodge duality operator on fc-forms in dimension n. 
Formula (| A.3|) also yields solutions on non-contractible domains. In that case, 
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one may allow singular stream- forms /. Equation (|A.3I) may also be replaced 
by the more general expression 

w = 1 i n i2\\ ( A - 4 ) 

for arbitrary closed (n — fc)-forms a. 

The Frobenius condition for the k forms ui expressed by (|A.3|) holds for 
1-forms r satisfying 

r A *df = d* df — d\ogp(^(\df\ 2 )) A *df . (A.5) 

This is equivalent to the existence of an integrating factor for uj when T can be 
made exact. This can always be done in the case in which w is a 1-form or an 
(n — l)-form. Equation (|A.5[) can be satisfied with 

r = r 1 -diog P (^(\df\ 2 )) , (A.6) 

whenever Ti satisfies 

d*df = T 1 A*df . (A.7) 

This argument shows that the ability to choose T exact does not depend on the 
particular type of density function employed, although the value of T itself docs 
depend on p. For ui expressed by (IA.4I) . formulas (IA.5I) , (|A.6[) and (|A.7[) hold 
with df replaced by a. 

The equations studied in Sees. [5] and [3] correspond to the cases of 1-forms 
in dimensions 2 and n, respectively. The system studied in Sec. [4] corresponds 
to the case of (n — l)-forms in dimension n. In these cases the differential 
forms can be easily interpreted as vector functions: a 1-form can always be 
interpreted as a vector function, as can an (n — l)-form, being the Hodge-dual 
of a 1-form. For the intermediate cases, in which k is neither 1 nor n — 1, the 
description in terms of differential forms is more natural. For this reason, an 
interesting - perhaps the most interesting — illustration of this Appendix is 
provided by simply taking k = 2, n = 4, and p as in (|1.9[) . This corresponds to 
the Born-Infcld model in dimension 4. Our method yields the solutions 

= * df : / = y a ( x ) dx J e A 1 (fi) . (A.8) 

The solutions uj+ are defined (and uniformly bounded) for smoothly prescribed 
generalized stream 1-forms /, while the solutions w_ require the additional 
condition \df\ > 1 and are unbounded for choices of generalized stream forms 
which satisfy \df\ — 1 at points of the domain f2. In this example, the Frobenius 
condition holds for 1-forms T satisfying 

ld(\df\ 2 ) , x 

rA * df = d * df -2\WTT- 
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